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Abstract 

We study the effect of theoretical and parametric uncertainties on the ability of future Higgs 
coupling measurements at the International Linear Collider (ILC) to reveal deviations from the 
Standard Model (SM). To quantify the impact of these uncertainties we plot Ax 2 = 25 contours 
for the deviations between the SM Higgs couplings and the light Higgs couplings in the m™ ax 
benchmark scenario of the Minimal Supersymmetric Standard Model (MSSM). We consider the 
theoretical uncertainties in the SM Higgs decay partial widths and production cross section and 
the parametric uncertainties in the bottom and charm masses and the strong coupling a s . We find 
that the impact of the theoretical and parametric uncertainties is moderate in the first phase of 
ILC data-taking (500 ftr 1 at 350 GeV centre-of- mass energy) , reducing the "reach" in the CP-odd 
MSSM Higgs mass M A by about 10% to ~ 500 GeV, while in the second phase (1000 ftr 1 at 1000 
GeV) these uncertainties are larger than the experimental uncertainties and reduce the reach in 
Ma by about a factor of two, from ~ 1200 down to ~ 600 GeV. The bulk of the effect comes from 
the parametric uncertainties in m& and a s , followed by the theoretical uncertainty in r^. 
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I. INTRODUCTION 



An important component of the physics case for the International Linear Collider (ILC) 
is its ability to perform high-precision measurements in the Higgs sector [l[ [2], [J. In addition 
to measuring the Higgs mass, its spin and CP quantum numbers, and its self-coupling which 
characterizes the shape of the Higgs potential, the ILC will measure the Higgs couplings to 
a variety of Standard Model (SM) particles. In the SM, the masses of these particles arise 
solely through their couplings to the Higgs, so that these couplings are predicted in terms of 
the corresponding SM particle masses. The SM Higgs production cross sections and decay 
branching ratios are then fixed once the Higgs mass has been measured. 

In extended models, however, the masses of the SM particles can receive contributions 
from more than one source, e.g., from couplings to each of two Higgs doublets. Further, the 
SM-like Higgs boson in extended models is typically an admixture of states from different 
scalar multiplets, leading to modifications of its couplings to SM particles by additional 
mixing angles. The high-precision ILC measurements of Higgs couplings thus provide a 
powerful tool for probing the structure of the Higgs sector and distinguishing the minimal 
SM from extended models. 

Previous studies have quantified the "reach" of the ILC measurements 

within the context of the Minimal Supersymmetric Standard Model (MSSM). In this paper 
we study the impact on this reach due to the theoretical and parametric uncertainties that 
enter the SM predictions for the Higgs cross section and decay branching ratios. A first 
attempt at including the parametric uncertainties from the strong coupling and the bottom 
and charm quark masses was made in Ref. |5|] by varying these parameters within their 
uncertainties and treating the resulting variation in the SM Higgs branching ratio predic- 
tions as an additional uncertainty to be added in quadrature with the uncertainty on the 
experimental measurement. An additional step was taken in Ref. {9] where the theoretical 
uncertainties in the SM 77 — > H — > bb rate were estimated based on the size of the computed 
radiative corrections. We refine and extend the analysis of the parametric uncertainties by 
taking into account their correlated effects on Higgs observables. We also include the full 
set of current theoretical uncertainties in the SM Higgs decay partial widths and production 
cross section, again taking into account their correlated effects on the ILC Higgs observables. 

With the theoretical and parametric uncertainties in hand, we must choose a measure to 
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quantify their impact. We choose to work in the m™ ax benchmark scenario of the MSSM as 



defined in Ref. 



10J. We scan over points in the M^-tan/3 plane and construct a chi-squared 



(x 2 ) observable between the SM and MSSM predictions. We plot contours of A% 2 = 25, 
corresponding to a 5cr discrepancy between the SM and the sample MSSM point. These 
contours allow us to show the impact of adding each of the theoretical and parametric 
uncertainties on top of the experimental uncertainties. Because the main purpose of this 
paper is to examine the impact of the theoretical uncertainties, we make no attempt to scan 
over more general sets of MSSM parameters or over non-supersymmetric models. 

We begin in the next section with a brief overview of the MSSM Higgs couplings, frame- 
work of our analysis. In Sec. Ill II we describe our \ 2 procedure for dealing with the uncorre- 
cted and correlated sources of uncertainty. In Sec. [IV] we collect the expected experimental 
uncertainties in the Higgs measurements and describe our treatment of the theoretical and 
parametric uncertainties, which typically feed in to multiple measurements in a correlated 
way. We provide a table of the dependence of individual Higgs decay partial widths on the 
input parameters m&, m c and a s so that our study can be updated in a straightforward 
way if the precision on these parameters improves. We present our numerical results in 
Sec. [V] Section [VI] is reserved for discussion and conclusions. The current state of the art 
for the radiative corrections to SM Higgs decays, Higgs production in e + e~ — > vi/H, and the 
bottom and charm quark mass extraction from low-energy experimental data is reviewed in 
Appendices El EJ and [Cj respectively. 



II. THEORETICAL FRAMEWORK 

The couplings of the lightest CP-even Higgs boson h° of the MSSM can be written at 
tree level in terms of the corresponding SM Higgs couplings as 



9h°tt 

9H SM it 


9h°cc 
9H S mcc 


= sin(/3 — 


a) + cot (3 cos(/3 — a) 


9h°bb 
9H SM bb 


9b? TT 

9H S mtt 


= sin(/3 — 


a) — tan j3 cos(/3 — a) 


9h°WW 


gh°zz 


= sin(/5 — 


a), 


9h sm ww 


9h sm zz 



(1) 



3 



where a is the mixing angle that diagonalizes the mass matrix for the CP-even states h° and 
H°, and tan/5 = i^/fi is the ratio of vacuum expectation values of the two MSSM Higgs 
doublets. At tree level we have 

cos(/3-a)~isin4/3^f, (2) 

where Ma is the mass of the CP-odd neutral MSSM Higgs boson. Note that in the decoupling 
limit Ma 3> rriz [llj], cos(/3 — a) goes to zero and the h° couplings approach their SM values. 
Using Eqs. (0Q) and (j2J) to expand the Higgs partial widths in powers of m z /M\, we obtain 
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cot/5sin4/3^|^-4cot 2 /3^|, (3) 

where the last equality in each line uses the large tan/5 approximation sin 4/3 ~ — 4cot/3. 
In particular, for large tan /3, r& and T r exhibit the largest deviations from their SM values, 
while T\y approaches its SM value very quickly with increasing Ma- 

Beyond tree level, the Higgs couplings receive radiative corrections, which mainly impact 
the CP-even Higgs mass matrix [12j . These can be taken into account by defining an effective 
mixing angle a e s for use in place of a in the above formulas. Significant radiative corrections 
can also appear in the h°bb vertex due to squark-gluino loops, with their dominant piece 



parameterized as A b 12J. Including this vertex correction, the h°bb coupling can be written 
as Q 

9hW r • ,o v , a (a ^ 1 - A fe cot a cot /3 

= [sm(/3 — a) — tan p cos(p — a)\ . (4) 

9H SM bb ' 1 + A b 

In the limit of large Ma, cot a cot (3 — > — 1 and the A& corrections also decouple. 

In this paper we make our x 2 comparisons using the implementations of the SM and 



MSSM Higgs couplings in the public Fortran code HDECAY 13J. HDECAY computes the h° 
couplings using the a e fr approximation, incorporating also the potentially significant A5 
corrections to the h°bb coupling. For the loop-induced couplings of h° to photon or gluon 
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pairs, HDECAY includes the shift in the contribution from SM particles in the loop 1 due 
to the modified couplings in Eq. (pQ), as well as the new contributions to the amplitude 
from supersymmetric particles running in the loop and from the charged Higgs boson for 
h° — ► 77. For comparisons involving MSSM Higgs production cross sections in WW fusion 
and ZH associated production, we scale the corresponding SM production cross section 
by 9h°vv/9H SM vv (V = W or Z) using the a e ff approximation of HDECAY. We neglect the 
additional non- universal SUSY electroweak radiative corrections (e.g., box diagrams) beyond 
this approximation. While a full MSSM analysis should include complete non-universal 
SUSY corrections to Higgs partial widths and production cross sections, the approximations 
used in HDECAY are sufficient for our purposes in this paper. 

III. CHI-SQUARED ANALYSIS 

To quantify the impact of the experimental and theoretical uncertainties on the ability 
of ILC measurements to reveal deviations from the SM, we construct a \ 2 according to 

n n 

* 2 = E E^ 1 - Q? a ) W(9f 1 - ( 5 ) 
i=i 3=1 

where Qi are the observables being compared between models M\ and M 2 , and [c 2 ]^ 1 is the 
inverse of the covariance matrix for the quantities Qi, defined according to 

m 

4 = SijUiUj + ^cfc*. (6) 

k=l 

Here Ui is the uncorrelated uncertainty in the observable Qi, represents the kth source 
of correlated uncertainty in Qi, and 5ij is the Kronecker delta. In the absence of correlated 
uncertainties, afj is a diagonal matrix with the (positive definite) squares of the uncorre- 
lated uncertainties Ui down the diagonal. Correlated uncertainties introduce off-diagonal 
terms, which can be positive or negative depending on the sense of the correlation between 
observables Qi and Qj. 

1 We added the charm quark loop to the SM T gg calculation in HDECAY, to be consistent with its inclusion 
in the MSSM calculation. 
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In our analysis the observables Qi consist of Higgs branching fractions and Higgs produc- 
tion cross sections times branching fractions (i.e., rates in a particular channel). We take 
the expected experimental uncertainties on these quantities at the ILC from the literature. 
Because the existing experimental studies have not quantified the correlations between the 
measured observables, we treat the experimental uncertainties as uncorrelated. The theo- 
retical and parametric uncertainties affect multiple observables, and we therefore treat them 
as correlated uncertainties. 



IV. INPUTS 

A. Experimental inputs 

We consider Higgs measurements from two stages of ILC running: 

• Phase 1: 500 fb _1 of integrated luminosity at 350 GeV centre-of-mass energy, with 
no beam polarization. Expected uncertainties in the SM Higgs branching ratios for 
Higgs masses m# = 120 and 140 GeV are taken from Ref. 14| and summarized in 
TableE 

• Phase 2: 1000 fb _1 of integrated luminosity at 1000 GeV centre-of-mass energy, 
with beam polarizations of —80% for electrons and +50% for positrons. Expected 
uncertainties in the SM Higgs production cross section times branching ratios for 
nin = 115, 120 and 140 GeV are taken from Ref. [15j and summarized in Table HT1 
For the Phase 2 analysis, we include the existing Phase 1 measurements in the x 2 in 
addition to the new measurements at 1000 GeV. 

The centre-of-mass energy of 350 GeV was chosen for the study in Ref. [l4j because it 
is near the peak of the e + e~ — > ZH cross section for the Higgs masses considered. It is 
also near the top quark pair production threshold, so that the Higgs data could be collected 
simultaneously with a top threshold scan. Running with polarized beams would improve 
the results by boosting the cross section and/or suppressing backgrounds. The actual ILC 
run plan will of course depend on what is discovered at the CERN Large Hadron Collider 
(LHC), and may include multiple threshold scans or running at the maximal first-phase 
design energy of 500 GeV. 
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SM Higgs branching ratio uncertainties from 500 fb 1 at 350 GeV 
m H = 120 GeV 140 GeV 



BR(66) 


2.4% 


2.6% 


BR(cc) 


8.3% 


19.0% 


BR(rr) 


5.0% 


8.0% 


BR(WW) 


5.1% 


2.5% 


BR{gg) 


5.5% 


14.0% 



TABLE I: Expected fractional experimental uncertainties in the SM Higgs branching ratios for 



mu = 120 and 140 GeV, from Ref. LJ] (500 fb 1 at 350 GeV centre-of-mass energy with no beam 
polarization). 



SM Higgs cross section times BR statistical uncertainties from 1000 fb 1 at 1000 GeV 
m H = 115 GeV 120 GeV 140 GeV 


a x BR(66) 


0.3% 


0.4% 


0.5% 


a x BR(WW) 


2.1% 


1.3% 


0.5% 


a x BR(gg) 


1.4% 


1.5% 


2.5% 


a x BR(77) 


5.3% 


5.1% 


5.9% 



TABLE II: Expected fractional experimental uncertainties in the SM Higgs production cross section 



times decay branching ratios for uih = 115, 120 and 140 GeV, from Ref. 15(] (1000 fb -1 at 1000 
GeV centre-of-mass energy). Uncertainties are statistical only. Beam polarizations of —80% for 
electrons and +50% for positrons are assumed. 



The study in Ref. [l5|] at 1000 GeV centre-of-mass energy found the expected statistical 
uncertainties on Higgs rates in various channels. The analysis selected events with e + e _ — > 
uuH, thus including Higgs production via WW fusion and via e + e~ — > ZH with Z — > vv. 
The Higgs decay products were then selected with the requirement that the visible energy 
in the event add up to the Higgs mass. The study in Ref. [lj]] evaluated the statistical 



uncertainties only. In our x 2 analysis we add an overall luminosity uncertainty of 0.1% [l] for 
the Higgs rate measurements at 1000 GeV, completely correlated among the four channels. 
In this high-energy phase of ILC running, Higgs production is dominated by WW fusion, 
the cross section for which grows with centre-of-mass energy for a fixed run- It is thus 
advantageous to run at the highest possible collider energy to maximize the statistics. Beam 
polarization was chosen to maximize the WW fusion cross section. 

In our analysis we scan over the M^-tan/? plane in the m™ ax scenario of the MSSM. At 
each parameter point we take the Higgs observables for h° and those for a SM Higgs boson of 
the same mass as the input for our x 2 - For Higgs masses between 120 and 140 GeV, we obtain 
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Theory uncertainty 
Higgs partial width in literature in HDECAY 





1% 


1% 




n mo/ 
U.Ul/o 


U.Ul /o 




0.5% 


5% 


r 


3% 


16% 


r 


0.1% 


4% 




4% 


4% 


Higgs production cross section 








0.5% 





TABLE III: Estimated fractional theoretical uncertainties in the SM Higgs partial widths and 
production cross section due to uncalculated higher order corrections. See Appendices \M and [B] 
for details. 

the expected experimental uncertainty on the Higgs observables using linear interpolation 
between the values given in Tables [I] and [TT1 In the m™ ax scenario the h° mass never exceeds 
130 GeV. For h° masses below 120 GeV, which occur in this scenario only at low tan f} values 
< 4, we use the experimental uncertainties for m# = 120 GeV from Tableland apply linear 
interpolation only to the values given in Table HT1 between m# = 115 and 120 GeV. 



B. Theoretical and parametric uncertainties 



1. Higgs decay partial widths 

The theoretical uncertainties in the SM Higgs decay partial widths due to uncalculated 
higher order radiative corrections are summarized in Table IIHI We estimated these numbers 
based on the size of the known higher order corrections for Higgs masses in our range of 
interest; for details see Appendix [S] In all cases we use the best theory uncertainty available 
in the literature for our x 2 calculation; in some cases the theory uncertainty in the HDECAY 
calculation is larger because not all available radiative corrections have been implemented 
into the HDECAY code Q. 

Working entirely in terms of fractional uncertainties, the correlated uncertainty cf in 
branching ratio i due to the theoretical uncertainty in the SM prediction for partial width 
k is given by the formula 
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Parameter Value Percent uncertainty Source 



a s {m z ) 0.1185 ±0.0020 1.7% [20] 
m- b (M b ) 4.20 ± 0.04 GeV 0.95% [21] 
rn c (M L ) 1.224 ± 0.057 GeV 4.7% [22] 

TABLE IV: Central values and uncertainties of a s and the MS bottom and charm quark masses 
rrTb(Mb) and m^(M c ) used in the analysis. 

where a^ k is the (fractional) theoretical uncertainty on partial width T^- The normalized 
derivatives are given analytically by 

T k dBR t ) -BR fc for % ± k (g) 



BRi dT k ] ft _ BRh ) for i = k 
2. Cross section for e + e~ — * vuH 

At 1000 GeV centre-of-mass energy the production cross section for e + e _ — > uuH is 
dominated by WW fusion, with a subleading contribution from e + e~ — > ZH with Z — > 
z/z/. The SM cross section is currently known including the one-loop electroweak radiative 



corrections 



13, 



18 



19| . We take the remaining theoretical uncertainty in the cross section 
to be 0.5%; see Appendix [B] for details. This cross section uncertainty is included in the 
X 2 calculation in the same way as the luminosity uncertainty discussed earlier, completely 
correlated among the four rate measurements. 



3. a s , m^Mb), and m c (M c ) 

The most important sources of parametric uncertainties in the SM Higgs coupling cal- 
culations are the strong coupling a s and the bottom and charm quark masses. The input 
values used in our analysis are summarized in Table IIVI For a s we use the current world 
average from the Particle Data Group [20| . For the bottom and charm quark masses we use 
the values from fits to the kinematic moments in inclusive semileptonic B meson decays; 
details are given in Appendix O We make the approximation that the MS quark mass 
m^(Mg) evaluated at the corresponding quark pole mass M q (used by HDECAY) is the same 
as the mass m^(m^) evaluated at its own running mass (extracted from the B decay fits). 

We implement the parametric uncertainties in the \ 2 calculation as follows. Again, all 
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uncertainties are fractional. The correlated uncertainty in branching ratio i due to the 
parametric uncertainty in Xj G {a s (mz), m&(Mb), m^(M c )} is given by 



Xj dBRi 
BRj dxj 



/c=i 



BRi dT k 



xj dT k 
T k dxj 



(9) 



where a Xj is the (fractional) uncertainty on the parameter Xj. The second equality was 
obtained using the chain rule; the first term inside the sum was given explicitly in Eq. (jSJ). 

The calculation thus consists of evaluating the derivatives dT k /dxj while holding constant 
all the parameters other than Xj. In principle, this can be done by varying each of the inputs 
Xj in HDECAY and reading off the resulting variation in IV This is complicated by the fact 
that HDECAY takes as inputs not the running masses ra&(Mft) and m^(M c ), but the pole 
masses Mb and M c themselves. Varying M& while holding M c constant is not the same as 
varying rn%(Mb) while holding m^(M c ) constant; likewise, varying a s (mz) while holding Mb 
and M c constant results in a variation of Wb(Mb) and m^(M c ). 

To deal with this, we use the chain rule and next-to-leading-order (NLO) approximations. 
We start by writing the quark pole masses as functions of our desired input parameters: 
Mb = Mb(mb, a s ) and M c = M c (Wb, rn~ c , a s ), where the running quark masses are understood 
to be evaluated at their corresponding pole mass scales and a s = a s (mz)- Any Higgs partial 
width T can then be expressed as T = T(M b (j?Tb,a s ), M c (mfc, m^, a s ), The chain rule 
then gives 



or 




or 




0M b 


drfib 


m c ,a 3 


0M b 


M c ,a a 


drrib' c 


dT 




or 




dM c 


dwr c 




dM c 


M b ,a 3 


dm~ c - 


or 




dT 




dT 


da s 


mij,m c 


da s 


M b ,Mc 


+ DM 



+ 



m b ,a s 



dT 



dM r 



dM c 

M b ,a s 



(10) 



m c ,a s 



dM h 



da. 



+ 



dT 



n>b 



dM r 



dM n 



M b ,a a 



da. 



m b ,m c 



HDECAY can be used in a straightforward way to numerically evaluate [dT/dMb]M c ,a s , 
[dT/dM c ]M b ,a a , and [dT/da s ]M b ,Mc- We did this by varying each of the input parameters 
Mb, M c , and a s about its central value over a range of approximately twice its uncertainty 
as given in Table HVl We also checked explicitly that the dependence of each of the partial 
widths was linear under these (small) variations. 
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We evaluate the remaining derivatives analytically as follows. The relation between the 
pole mass M q and running mass = r?T q ~(M q ) at NLO is 



M„ = m n 



1 + —a s (M q ) 



Differentiating this with respect to m q yields 



dM q _ . L _£ 4 da s (M q ) 

dmZ 3tt s 9 9 3ir dmZ 



(11) 



(12) 



The derivative of a, can be found from the NLO formula 



a s (M q 



a s (n) 



1 + (6/27r)a a (/x) log(M 9 //i) 



(13) 



where b = 11 — 2Nf/3 and Nf is the number of active quark flavours in the running of a s 
between the scales /x and M q . In particular, da s (M q )/drrTg~ = 0(a 2 s ); working to NLO we 
will neglect this term. The derivative of interest then becomes 



dM q 1 4 M„ 
om q 6 m q 



(14) 



i.e., the normalized derivative (m q / M q ) [dM q j dm q ) = 1 at NLO. A similar calculation yields 
dM c /dm = at NLO. 

The derivatives of the quark pole masses with respect to a s are evaluated as follows. 
From Eq. (fTTj) we have, to NLO, 



_ 4 da s {M q 



da s (m z ) 



m r , 



37r da s (mz) 



(15) 



Equation (|T3|) can be used to evaluate the a s derivative at NLO, 



da s (M q ) 
da s (m z ) 



a s {M q ) 
a s {m z ) 



(16) 



2 This formula automatically takes into account the change in the number of quark flavours in the running 
at the b threshold. 
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Normalized derivatives of Higgs partial widths 





f~\l ( f~H rr 1 


rnZ( AAA 


wr( M \ 




120 GeV 140 GeV 160 GeV 


120 GeV 140 GeV 160 GeV 


120 GeV 140 GeV 160 GeV 


^bb 


-1.177 -1.217 -1.249 


2.565 2.567 2.568 


0.000 0.000 0.000 


r - 

1 cc 


-4.361 -4.400 -4.432 


-0.083 -0.084 -0.084 


3.191 3.192 3.192 


1 99 


2.277 2.221 2.175 


-0.114 -0.112 -0.104 


-0.039 -0.032 -0.027 


v 

i 77 


0.002 0.002 0.001 


0.010 0.008 0.005 


0.012 0.009 0.005 



TABLE V: Normalized derivatives, (x/T)(dT/dx), of SM Higgs partial widths with respect to the 
parameters a s (mz), m&(Mfe), and m^(M c ) for m# = 120, 140, and 160 GeV. See text for details. 
The central values of the three input parameters were taken from Table IIVI 

Combining Eqs. (I15p and (I16p we find for the normalized derivative, 



aAmz 



M q da s (m z ) 



M q 3n a s {m z ) 



(17) 



where in the last step we keep only terms up to NLO. 

The dependences of the SM Higgs partial widths Tj on variations in the inputs Xj G 
{a s , frTb(Mb), 7 m c ~(M c )} are summarized in Table IVl in the form of normalized derivatives, 
(xj/Ti)(dTi/dxj), evaluated at the central values of the inputs Xj. For example, for m# = 
120 GeV, a 1% increase in Wib(Mb) leads to a 2.6% increase in T bb , while a 1% increase in 
OL s {mz) leads to a 1.2% decrease in T bb : The dependence of the normalized derivatives on 
the Higgs mass is not very strong. We again use linear interpolation to find the appropriate 
normalized derivatives for Higgs mass values between those given in Table |V] Finally, we 
note that the Higgs partial widths to WW, ZZ, and lepton pairs, and the production cross 
section for e + e~ — > vvH, do not depend on a s , Tnb(M b ), or m^(M c ). 

The general pattern of the normalized derivatives in Table IVl can be understood as follows. 
The H —* qq partial widths are proportional to rn q (rriH) [see Eq. flAlp ] . Neglecting QCD 
running of the quark mass, the normalized derivative of T qq with respect to should thus 
be equal to 2. The values in Table |V] are actually somewhat larger than 2 because of the 
effect of QCD running, which causes m^(n) to decrease as the scale /i increases. Raising 
rn q ~(M q ) by a small amount reduces the range of scales (from M q to over which the QCD 
running is applied, so that m^(m//) is even larger than it would be just from the increase in 
rn q ~(M q ) without the effect of the running. This effect is more pronounced in T c5 than in T bb 
because a s is larger at the charm quark mass scale than at the bottom mass scale, resulting 
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in faster running for rn~ c near the scale M c . Similarly, increasing a s (mz) strengthens the 
QCD running of W q , resulting in a reduction of Y qq . Changing a s also affects the QCD 
corrections to T qq ; however, the dominant effect is absorbed already into the running quark 
mass. 

The slight dependence of r cg on nTb(Mi,) is due to the effect of the b threshold on the 
running of rrT c . Adding the b quark to the renormalization group equations slows down the 
running of rn~ Cl so that increasing the b mass slightly reduces the high-scale value of WT C and 
thus reduces T c5 . Note that by working to NLO in evaluating the derivatives above, we 
have neglected contributions to the nib dependence of T c5 of the same order as those that 
give the dependence shown in Table |V] a higher order treatment would thus give different 
results. However, the numerical effect of these terms is negligible, so we do not attempt a 
more precise evaluation here. 

Y gg is proportional to a 2 s at leading order. NLO QCD corrections increase Y gg rather 
significantly, leading to a normalized derivative somewhat larger than 2. T gg is dominated 
in the SM by the top quark loop, so that the effect of varying rfTb(Mb) and m^(M c ) is small. 
The anticorrelation in this case is due to destructive interference between the real parts of 
the b and c quark loop amplitudes and the top quark loop amplitude. 

Finally, the similar sensitivity of T 77 to m^M;,) and m^(M c ) can be understood by noting 
that the bottom and charm quark loops contribute mainly through their interference with 
the dominant top and W boson loops, and that their contributions go like mb(mH)Ql and 
m^(m H )Ql, respectively, while m^(fi) /m^(fi) ~ 0.235 ±0.012 [2l|] and Q 2 C = AQl so that the 
two contributions are actually comparable. 

V. RESULTS 

Our main results are summarized in Fig. [H where we plot A% 2 = 25 contours for the 
deviations between the SM Higgs and the lighter CP-even Higgs in the m™ ax scenario of the 
MSSM. To the left of these contours, the MSSM Higgs can be distinguished from the SM 
Higgs at more than the 5a level. The dashed curves show the reach including experimen- 
tal uncertainties only, while the solid curves show the effect of adding the theoretical and 
parametric uncertainties. 

For the first phase of ILC running (the left pair of curves in Fig. [TJ 500 fb _1 at 350 
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400 600 800 1000 1200 1400 



M A [GeV] 

FIG. 1: Contours of Ax 2 = 25 from the experimental uncertainties only (dashed lines) and includ- 
ing all theoretical and parametric uncertainties (solid lines). The pair of contours on the left are 
for Phase 1 and the pair on the right are for Phase 2. 

GeV, experimental uncertainties from Tabled]), the impact of the theoretical and parametric 
uncertainties is to reduce the reach in Ma by about 50 GeV, or roughly 10%. Details 
are shown in Figs. [2] and [21 where we plot separately the contribution of each source of 
parametric and theoretical uncertainty. The bulk of the effect is due to the parametric 
uncertainties in mb(Mf,) and a s (mz), as shown in Fig. [2j Smaller contributions come from 
the theoretical uncertainty in r&, and for low tan/5 < 6 from the parametric uncertainty in 
m^(M c ) and the theoretical uncertainty in T g . The change with tan/3 is partly due to the 
decrease in rrih with decreasing tan/3, as shown by the horizontal dotted lines in Figs. [2] and 

El 

The relative importance of the different branching fraction measurements in the first phase 
of ILC running is illustrated in Table IVIl We list the fractional deviations of the MSSM 
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M A [GeV] 

FIG. 2: Contributions of individual sources of parametric uncertainty to the A% 2 = 25 contours 
in Phase 1. 





Phase 1 sample point: Ma = 537.6 GeV, tan/3 


= 20 




Observable 


Shift 


Expt uncert 


Pull 


Thy+param uncert Total uncert 


Pull 


BR(66) 


8.1% 


2.5% 


3.25 


1.6% 


3.0% 


2.71 


BR(cc) 


-12.0% 


13.2% 


-0.90 


16.1% 


20.8% 


-0.57 


BR(rr) 


10.0% 


6.4% 


1.56 


1.8% 


6.6% 


1.51 


BR(WW) 


-11.6% 


3.9% 


-2.96 


1.8% 


4.3% 


-2.68 


BR{gg) 


-14.7% 


9.4% 


-1.56 


5.8% 


11.1% 


-1.33 



TABLE VI: Contributions to the A% 2 at a Phase 1 sample point lying on the A% 2 = 25 contour 
including only experimental uncertainties in Figs. [2] and [3j Shown are the fractional deviations in 
Higgs branching ratios in the m™ ax scenario compared to the SM, the experimental uncertainties 
for the corresponding Higgs mass = 129.2 GeV taken from a linear interpolation in Higgs mass 
between the values in Table [U the resulting "pulls," i.e., the branching ratio deviations normalized 
by their experimental uncertainties, the combined theoretical and parametric uncertainties on 
the observables, the total uncertainty, and the resulting pull from the total uncertainty. Not 
shown are the off-diagonal elements of the correlation matrix due to the theoretical and parametric 
uncertainties. 
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FIG. 3: Contributions of individual sources of theoretical uncertainty to the A% 2 = 25 contours in 
Phase 1. The additional theory uncertainties not plotted separately have a negligible effect. 

Higgs branching ratios from their SM values at the "Phase 1 sample point" Ma = 537.6 GeV 
and tan/? = 20, which lies on the A% 2 = 25 contour from experimental uncertainties only 
and yields = 129.2 GeV. The relevant branching fractions deviate by roughly 10% from 
their SM values at this sample point. The biggest contributions to the A% 2 come from the 
best-measured branching fractions, BR(66) and BRfWW 7 ). Including the theoretical and 
parametric uncertainties moderately degrades the precision in all the channels, reducing the 
Ax 2 from 25 down to 17.4 at the sample point. This A% 2 is slightly worse than the value 
of 18.9 that would be obtained by summing the squares of the pulls in Table [VU this is 
due to the effect of the off-diagonal elements in the correlation matrix. The effects of the 
most important parametric uncertainties, in m^M^) and a s (rriz), are anticorrelated between 
BR(66) and BR(W / W / ) and thus can mimic the anticorrelated shift in these observables that 
occurs in the MSSM. 
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FIG. 4: Contributions of individual sources of parametric uncertainty to the A% 2 
in Phase 2. 



25 contours 



For the second phase of ILC running (the right set of curves in Fig. [TJ combining mea- 
surements from 1000 fb _1 at 1000 GeV and 500 fb" 1 at 350 GeV, experimental uncertainties 
from Tables HI and HT1) . the experimental uncertainties are considerably smaller. The theo- 
retical and parametric uncertainties thus have a much more significant impact: their effect 
is to reduce the reach in Ma by about a factor of two, from about 1200-1400 GeV down to 
about 525-700 GeV. Details are shown in Figs. H] and [5j The dominant contributions by far 
are the parametric uncertainties in m^M^) and a s (mz) (Fig. H]), followed by the theoretical 
uncertainty in Tf, and for low tan (3 < 4 in T g (Fig. [5]). 

The impact of the parametric and theoretical uncertainties on the individual channels 
used in the Phase 2 A% 2 calculation is illustrated in Table IVHI We list the fractional devi- 
ations of the MSSM Higgs observables from their SM values at the "Phase 2 sample point" 
Ma = 1302.4 GeV and tan/3 = 20, which lies on the A% 2 = 25 contour from experimental 
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FIG. 5: Contributions of individual sources of theoretical uncertainty to the A% 2 = 25 contours in 
Phase 2. The additional theory uncertainties not plotted separately have a negligible effect. 



uncertainties only and yields rrih = 129.3 GeV. The relevant observables deviate by roughly 
2% to 4% from their SM values at this sample point, except for crxBR(77) which is very 
close to its SM value. Ignoring the parametric and theoretical uncertainties, the biggest 
contribution to the A% 2 comes from the best-measured observable, a x BR(66), followed by 
a x BK(gg) and o x BRfWVT). The contribution of the branching fraction measurements 
from 350 GeV running to the A% 2 is small. 

Including the theoretical and parametric uncertainties severely degrades the precision of 
the three most important measurements, <rxBR(66), crxBR(W / W / ), and axBK(gg), reducing 
the Ax 2 dramatically from 25 down to 1.7 at the sample point. This A% 2 is somewhat 
worse than the value of 3.2 that would be obtained by summing the squares of the pulls 
in Table IVII( this is again due to the effect of the off-diagonal elements in the correlation 
matrix, including the correlations between the Phase 1 branching ratio measurements and 
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Phase 2 sample point 




= 1302.4 GeV, tan/3 = 


20 




Observable 


Shift 


Expt uncert 


r Ull 


Thy+param uncert Total uncert 


r Ull 


BR(bb) 


1.7% 


2.5% 


U.O I 


1.7% 


3.0% 


U.00 


BR(cc) 


-2.5% 


13.3% 


— U.ly 


16.1% 


20.8% 


n 1 o 
— U.lz 


BR(rr) 


2.1% 


6.4% 


U.o4 


i qo/ 

1.8/0 


R «07 

6.670 


n qo 




O 1 0/ 

— z.l/o 


q no/ 
o.y/o 


— U.Oo 


1.8% 


4.3% 


— U.4o 


BR(^) 


-4.6% 


9.4% 


—0.48 


5.8% 


11.1% 


—0.41 


a x BR(bb) 


1.7% 


0.45% 


3.72 


1.7% 


1.8% 


0.93 


j x br(ww) 


-2.1% 


0.93% 


-2.22 


1.9% 


2.1% 


-0.98 


a x BR(##) 


-4.6% 


2.0% 


-2.32 


5.8% 


6.2% 


-0.74 


a x BR(77) 


0.27% 


5.5% 


0.05 


1.9% 


5.8% 


0.05 



TABLE VII: Contributions to the A% 2 at a Phase 2 sample point lying on the A% 2 = 25 contour 
including only experimental uncertainties in Figs. [4] and 03 Shown are the fractional deviations 
in Higgs observables in the m™ ax scenario compared to the SM, the experimental uncertainties 
for the corresponding Higgs mass m/j = 129.3 GeV taken from a linear interpolation in Higgs 
mass between the values in Tables U and HU the resulting pulls, the combined theoretical and 
parametric uncertainties on the observables, the total uncertainty, and the resulting pull from the 
total uncertainty. Not shown are the off-diagonal elements of the correlation matrix due to the 
theoretical and parametric uncertainties. 

the Phase 2 rate measurements. 

VI. DISCUSSION AND CONCLUSIONS 

We studied the impact of the theoretical and parametric uncertainties in the SM predic- 
tions for Higgs production and decay on the ability of the ILC to reveal deviations from 
the SM Higgs. To quantify the impact of the theoretical and parametric uncertainties, we 
compared the SM predictions for Higgs observables to those in the m™ ax benchmark scenario 
of the MSSM. We plotted A% 2 = 25 contours for the deviations between the SM and MSSM 
Higgs observables both with and without the theory uncertainties, given the expected pre- 
cisions on Higgs measurements at a future ILC. We found that the impact of the theoretical 
and parametric uncertainties is moderate in the first phase of ILC data-taking (500 fb _1 
at 350 GeV centre-of-mass energy), reducing the reach in the CP-odd MSSM Higgs mass 
M A by about 10% to ~ 500 GeV, while in the second phase (1000 fb" 1 at 1000 GeV) these 
uncertainties are larger than the experimental uncertainties and reduce the reach in Ma 
by about a factor of two, from ~ 1200 down to ~ 600 GeV. The bulk of the effect comes 
from the parametric uncertainties in mj and a s , followed by the theoretical uncertainty in 
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IV The theoretical uncertainty in T g is also important for lower Higgs masses below about 
120-125 GeV. 

The single most important source of theoretical or parametric uncertainty is the bottom 
quark mass. In our analysis we used the measurement of rrib from a global fit to inclusive 
semileptonic B meson decay spectra 21J (see Appendix O for details). This is currently the 
best measurement of m&; its uncertainty of 0.95% is dominated by the theoretical uncer- 
tainties in the fit. The large impact of this parametric uncertainty on Higgs observables is 
due to (i) the effective 2.6 power with which rrTb(Mb) enters r&, (ii) the importance of the bb 
final state in ILC Higgs measurements, and (in) the fact that the deviations of the MSSM 
Higgs couplings from their SM values show up largely in the b and r sectors. The impact of 
the uncertainty in rrib on the long-term ILC Higgs program highlights the need for further 
theoretical work on the rrib extraction from B meson observables. 

The second most important source of theoretical or parametric uncertainty is a s (mz)- 



In our analysis we used the world-average PDG value [20| with an uncertainty of 0.0020 or 
1.7%. This measurement is expected to be improved by at least a factor of two at the ILC 
through precision measurements of event shape observables, the cross section ratio a t i/ 0^+^- 
above the ti threshold, and r^ ad /r^ pt at the Z pole (via the GigaZ option) [l|. We illustrate 
the effect of improving the uncertainty on a s (mz) to 0.0009 or 0.76% l| in Fig. [6j By itself 
this improvement has a relatively small effect because the uncertainty is still dominated by 
rrib, however, if the rrib extraction becomes significantly more precise, the improvement in 
a s will be valuable. 

The measurement of BR(cc) did not contribute significantly to the A% 2 because of its rel- 
atively poor experimental precision, particularly at the higher Higgs masses ~ 125-130 GeV 
that appear over most of the m™ ax parameter space. Even if the ILC measurement were 
to be improved, however, the parametric uncertainties in m^(M c ) and a s (mz) lead to a 
large uncertainty in the SM prediction for BR(cc) - 16% at our sample points - limiting the 
usefulness of this mode for model discrimination. The charm mass is by far the dominant 
uncertainty here, contributing 14.6%, while a s contributes 6.4%. Note however that the 
theoretical improvements in B meson decays required to reduce the uncertainty in rrib will 
also reduce the uncertainty in m c . 

We also note here that the theoretical calculation of higher-order QCD corrections to 
H — ► gg includes processes in which one of the final-state gluons splits into a cc or bb pair. 
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FIG. 6: Effect on the Phase 2 Ax 2 = 25 contours of improving the uncertainty on a s {mz) from 
1.7% to 0.76%, as projected for GigaZ in Ref. [J. 

Determining which part of these processes should be included in the H — > gg branching 
fraction and which in the H — ► cc or bb branching fractions will require interaction between 
the theoretical and experimental studies. The impact of this splitting on the bb final state 
is small, but it could potentially shift the numbers of events tagged as cc or gg by tens of 
percent (see the end of Appendix IA 41 for details). We have not attempted to address this 
question here. 

To understand the effect of the correlated uncertainties, we recall that for large tan /3, 
Tb and T T exhibit the largest deviations from their SM values, while Ty^ approaches its SM 
value very quickly with increasing Ma [Eq. (j3J)]. The largest parametric and theoretical 
uncertainties appear in the hadronic decay widths of the Higgs. At the same time, the most 
precisely measured observables in Phase 2 of the ILC are a x BR(fefe) and a x BR(gg), which 
suffer directly from the parametric and theoretical uncertainties, and a x BR(W / W / ), which 
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is affected indirectly through the Higgs total width. Adding a high-precision measurement 
of a x BR(rr) would allow a tight constraint on the MSSM Higgs unpolluted by the dom- 
inant uncertainties. This mode is affected by the parametric and theoretical uncertainties 
indirectly through the Higgs total width in the same way as o x BR(W / W / ), so that the ratio 
of these two non-hadronic modes is relatively clean. Further, T T exhibits a large deviation 
from its SM value in the MSSM [see Eq. (J3J)], leading to good potential sensitivity to the 
MSSM nature of the Higgs. 

To illustrate the potential improvement, we first note that for m# between 120 and 
140 GeV, the experimental uncertainties in the bb, WW, and gg final states improve by 
factors of 5-6, 4-5, and 3.5-5.5, respectively, between Phase 1 and Phase 2 (see Tables [I] 
and [III)- We thus expect that an improvement in the rr final state in Phase 2 by a factor 
of about four should be reasonable, and show the impact of a measurement of a x BR(rr) 
with an uncertainty of 1.3% (2.0%) at m H = 120 (140) GeV in Fig. El The effect of such 
a measurement on the reach in Ma including only experimental uncertainties is minor, 
because the precision assumed here is not particularly high compared to that of o x BR(66) 
and a x BR(H / H / ). However, the impact on the reach including parametric and theoretical 
uncertainties is quite significant, improving the reach from ~ 600 GeV to ~ 800 GeV. A 
measurement of o x BR(rr) would require a different selection than used in Ref. [15j], which 
required the visible energy in the event to add up to the Higgs mass; in the case of H — ► rr, 
some energy is lost through the neutrinos in the r decays. 

As the experimental studies for the ILC Higgs measurements are refined, we hope that the 
correlations among the experimental uncertainties in different channels will be quantified. 
In particular, we expect the hadronic decay modes H — > gg, bb and cc to be correlated, 
because the separation of these three channels relies on bottom and charm tagging. Because 
the shifts in Higgs observables in particular models are correlated due to the model structure, 
we expect the correlations in the experimental uncertainties to be important for evaluating 
model distinguishability. Once these correlated uncertainties are available, including them 
in a x 2 framework will be straightforward. 

We used the m™ ax benchmark scenario of the MSSM to quantify the impact of the theo- 
retical and parametric uncertainties. Our analysis can be extended in a straightforward way 
to Higgs coupling comparisons in other scenarios of the MSSM and to other beyond-the-SM 
Higgs sectors. Such a study would be interesting because the characteristic deviations from 
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FIG. 7: Effect on the Phase 2 Ax 2 = 25 contours of adding a measurement of a x BR(rr) with a 
precision of 1.3% (2.0%) at m H = 120 (140) GeV. 

the SM in different extended models may appear in different channels than for the MSSM 
m™ ax scenario, so that the relative importance of the various Higgs measurements will differ. 
Further, the shifts in the observables predicted in different extended Higgs models may pull 
in different directions relative to the correlations in the uncertainties in the SM predictions. 
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APPENDIX A: THEORETICAL UNCERTAINTIES IN HIGGS DECAY PAR- 
TIAL WIDTHS 

The radiative corrections to Higgs decays have been reviewed in Ref. [if]]; we give here 
a brief sketch of the known corrections. The theory uncertainties used in our analysis were 
summarized in Table HTT1 



1. H - qq 



The QCD corrections to H — > qq are known up to three loops (N 3 LO) 23, |24 L [25j for 
massless final-state quarks. A compact formula in the MS scheme is given in Ref. 



J25|| 
3: 



T(H - qq) 



3GVm# • 



AV2n q 



mJm H ) [Aqcd + A 



t > 



(Al) 



where 



Aqcd = 1 + 5.67^1 + (35.94 - 1.36iV F ) 

7T \ 7T / 



+ (164.14 - 25.77iV> + 0.259iV|) 



,(m H ) V 

7T J 



A,. 



1 + 0.20 + 0.04 + 0.003 
a s [m H - ' 1 r 



7T 



r _ 2 mjj 1 2 m K 
1-57 - - log + - log ^-^5 — 



3 ""° M 4 2 9 



0.02. 



(A2) 



In this calculation, the MS running quark masses fn^(M q ) evaluated at the quark pole mass 
M q are taken as a starting point and then run up to the Higgs mass scale using the standard 
QCD running in order to absorb large logarithms. The coefficient functions for the running 
of the quark masses are known to four- loop accuracy 26(|, as is the beta-function for the 
running of a s 271 ] . We thus estimate the uncertainty in the partial width from the running 
of m q (rriH) to be well below the 0.01% level. 

The perturbative expansion for Aqcd i n Eq- (IA2p is well under control; A t appears at 
NNLO and is roughly the same size as the NNLO term in Aqcd- We thus estimate the 
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uncertainty from uncalculated higher-order QCD corrections to be smaller than the last 
calculated term in each series: ~ 0.1% from Aqcd and ~ 1% from A t . 

We note here that an additional contribution to H — > qq can come from H — > gg* with 
the off-shell gluon splitting to qq; we have not included these effects in our calculations. See 
Sec. IA4I for further discussion. 

For Higgs decays to heavy quarks near the threshold region, the quark pole mass depen- 
dence must be included through the usual 3 = (1 — AM 2 /m 2 H ) 3 / 2 factor. The NLO QCD 
'or massive quarks are known 23|, as is the 0(Nf<y 2 ) piece of the massive NNLO 



corrections 
corrections 
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. HDECAY uses a linear interpolation between the NLO massive and the N 3 LO 
massless results near the threshold region [16]. For m# ~ 120 GeV, the quark mass effects 
decrease the H —>■ bb (cc) partial width by about 0.7% (0.1%). We estimate the remaining 
uncertainty from quark mass effects in the QCD corrections to be below the 0.01% level. 

The electroweak corrections to the Higgs decay into fermion pairs are known at one 
loop 29|, |20j; a convenient approximation formula is given in Refs. 16|, [311] • The QED 
corrections include a large logarithm log(m 2 H /M 2 ), which can be absorbed into the running 
fermion mass as in QCD. The electroweak corrections also include a top-mass-dependent 
piece which increases the H —>■ qq partial width by about 2.2% for decays to cc and by 
0.3% for decays to bb (the difference is due to the vertex correction involving a top quark 
in the case of bb, which partly cancels the universal correction). In addition, the two- and 
three-loop QCD corrections to the top-mass-dependent piece of the electroweak correction 



are known 
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331 ]. Their numerical size depends on the scheme chosen to define the running 



a s and top quark mass 33(; a suitable choice gives a relative correction to the top- mass- 
dependent piece of the NLO electroweak correction of 8% (two-loop) and 1% (three-loop). 
The remaining electroweak corrections are comparable in size, i.e., at the percent level. 
We thus estimate the remaining theoretical uncertainty from higher-order electroweak and 
mixed electroweak/QCD corrections to be smaller by a loop factor ~ a/ft, i.e., at the level 
of 0.01%. 

For large Higgs masses, electroweak corrections due to the Higgs self-coupling become 



relevant. These were calculated in Ref. 



34( and can increase T qq - by about 2% for m# ~ 1 



TeV. For a 120 GeV Higgs, they increase T qq - by about 0.2% at NLO and by another —0.002% 
at NNLO. Any remaining uncertainty from higher orders in the Higgs self-coupling is thus 
totally negligible in the intermediate Higgs mass range. 
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We thus estimate the theory uncertainty in T b g and T c5 to be at the 1% level, primarily 
from the At piece of the QCD corrections. 



2. H 



£+t 



Higgs decays to leptons receive the same electroweak corrections as the decays to quarks 
discussed in the previous section. As mentioned previously, we estimate that the theoretical 
uncertainty is at the level of two-loop electroweak corrections ~ (a/rc) 2 ~ 0.01%. 



3. H -> WW, ZZ 

For Higgs masses below the WW threshold, decays with one or both gauge bosons off- 
shell are important, with branching ratios for the doubly off-shell decays reaching the percent 
level for Higgs masses above about 100-110 GeV 16|. Until recently, calculations for off-shell 
vector bosons were only available at tree level, and radiative corrections to Higgs decays to 



WW and ZZ were known only in the narrow-widt 



i approximation. In this approximation, 



the one-loop electroweak corrections [29|, [35|, |36|, |37f amount to about 5% or less in the 

n 

intermediate mass range [16j. Furthermore, the QCD corrections to the leading top-mass- 



enhanced electroweak corrections of 0{a s Gpm 2 ) have been calculated at two |38l . |39| and 



three 



3, 



41] loops; these mixed electroweak/ QCD corrections amount to roughly 0.5%. 



HDECAY takes into account both singly off-shell and doubly off-shell decays at tree level but 
neglects these electroweak and mixed corrections [161 ]. 

Recently the complete NLO strong and electroweak radiative corrections for H 
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43] 



WW/ZZ — > 4/ have been calculated including the effects of off-shell gauge bosons 
The calculation was improved by including final-state photon radiation off the fermions, re- 
summed using a structure function approach, as well as higher-order corrections due to the 



461 ]. which increase the decay width 



Higgs self-coupling relevant for large Higgs masses [451 . 
for a 120 GeV Higgs by about 0.2% at NLO and by another 0.003% at NNLO. The complete 
calculation has been implemented in the Monte Carlo generator PR0PHECY4F 42j, |43j . For 
Higgs masses below 400 GeV, the theoretical uncertainty is below 1%; for our study we 
estimate the theory uncertainty in Ty/w and Yzz in the intermediate Higgs mass range to 



be at the 0.5% level 



44]. 
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4. H^gg 

The decay H — > gg appears first at one-loop mediated by a quark triangle. The top 
quark loop dominates the amplitude because it contributes proportional to the large top 
Yukawa coupling y t = gm t / \/2m w ~ 1. The bottom quark loop is also included in HDECAY; 
it reduces the H — > gg partial width by about 11% for m# = 120 GeV compared to including 
the top loop alone. We note that HDECAY v.3.103 [13] (downloaded June 2006) includes the 
charm quark loop in the calculation of the MSSM h° — > gg partial width, but not in the 
calculation of the SM H — > gg partial width. For consistency, we add the charm loop into 



the SM H —> gg calculation; it reduces Y gg by about 2% for m# = 120 GeV 



3 



The QCD corrections to H — > gg are most easily calculated in the heavy quark limit, 
in which m t — > oo is taken in the leading-order top triangle diagram (the amplitude goes 
to a constant in this limit). This allows the top triangle to be shrunk to a point and 
replaced with an effective coupling of the Higgs to two gluons. HDECAY includes the NLO 
QCD correction to T gg in this heavy quark limit only; in the intermediate Higgs mass range 
this correction shifts T gg upward by about 60-65% relative to the leading-order (LO) partial 



width. TheN 



O QCD corrections with full dependence on the mass of the quark in the loop 



are known |47|; in the intermediate Higgs mass range, including the full mass dependence 



shifts T gg upward by an additional 5% of the LO width 161 ] . This mass- dependent part of 
the correction is not included in HDECAY. 



The NNLO QCD correction is known only in the heavy quark limit 48]; it adds an 
additional 20% to the LO result. Very recently, the N 3 LO QCD corrections to T gg have 
been computed [49(], again in the heavy quark limit; they add another 2% to the LO result. 
By varying the renormalization scale between m^/2 and 2m#, a measure of the theoretical 
uncertainty of the calculation can be estimated: the scale dependence decreases from 24% 



(LO) to 22% (NLO) to 10% (NNLO), and finally to only 3% at N 3 LO [49J. Neglecting the 
mass dependence of the NLO correction and neglecting the NNLO and N 3 LO corrections 
entirely leaves HDECAY with about a 16% relative uncertainty on the calculated NLO partial 
width. 



The electroweak corrections of 0{GpmJ) to H — > gg were calculated in 50j in the heavy 



This number is for NF-GG = 5, i.e., the gluon splitting to bb or cc at NLO is not subtracted off T gg and 
added on to T b i cg ; see the end of this subsection for further discussion of this issue. 
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top approximation and they give a simple multiplicative factor which increases the gluonic 
decay width by about 0.3%. We thus expect the NNLO electroweak corrections and the NLO 
finite-top-mass effects to be below the per-mille level. This NLO electroweak correction is 
neglected in HDECAY. 

We thus take the remaining theory uncertainty in T gg to be roughly 3%, corresponding 
to the scale uncertainty in the N 3 LO QCD calculation. 

Finally, we note here a potential issue for the theoretical interpretation of Higgs partial 
width measurements that requires further study. Theoretical calculations of the higher-order 
QCD corrections to H — > gg include processes in which one of the final-state gluons splits 
into a quark-antiquark pair. The HDECAY code contains a switch that allows the user to 
specify the number of final-state quark flavours included in T gg ; setting NF-GG = 5 includes 
the five light quark flavours in T gg , while NF-GG = 4 separates off the contribution from 
gluon splitting to bb and adds it in to r 6 g. NF-GG = 3 does the same for gluon splitting to 
cc as well. 

Compared to NF-GG = 5 and for ran = 120 GeV, subtracting off the gluon splitting to 
bb and adding it to T b i (i.e., NF-GG = 4) raises T b i by about 1% and reduces T gg by about 
9%. Subtracting off the gluon splitting to cc and adding it to r cg (i.e., NF-GG = 3) raises 
r c g by almost 30% and reduces T gg by a further 12%. These effects are quite significant 
compared to the expected experimental precision in the gg and cc channels. The amount of 
this gluon splitting correction that should be subtracted off of T gg and included in r 6 g or T c5 
is an experimental question requiring a more detailed simulation of flavour tagging of Higgs 
decays into three-body final states. In our study we set NF-GG = 5 throughout and ignore 
this issue. 



5. H — ► 77 

The SM H — > 77 decay partial width receives QCD corrections to the one-loop diagrams 
involving quarks. Because the external particles in the ^77 vertex are colour neutral, the 

virtual QCD corrections are finite by themselves. The QCD corrections to the top quark 

[ I 

loop in T 77 are known analytically at NLO [51[ and as a power expansion up to third order in 



run/rut at NNLO 521 ] . The NLO corrections are of order 2% for m# < 2m w, and the NNLO 



corrections are at the level of a few per mille. These QCD corrections are neglected 16j in 
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HDECAY. 



T 77 also receives electroweak corrections, the complete set of which have been calculated 
at NLO neglecting the Yukawa couplings of all fermions except the top quark. The complete 
calculation consists of the two-loop diagrams containing light fermion loops and W or Z 
bosons (with the Higgs coupled to the W or Z, because the light fermion Yukawa couplings 
are neglected) 52 ]; the two-loop diagrams involving third-generation quarks and electroweak 



bosons 



3, 
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561 ]; and the purely bosonic corrections [56|]. Together, these NLO electroweak 
corrections are between —4% and 0% for 100 GeV < m# < 150 GeV. These electroweak 
corrections are neglected in HDECAY. 

The corrections to the part of the H — > 77 amplitude involving the Higgs coupling to 
light quarks (e.g., b) have not been computed. However, since the leading order contribution 
of the b quark loop to T 77 is at the few percent level, and the QCD and electroweak radiative 
corrections are themselves at the few percent level of the leading order piece, we estimate 
these uncalculated two-loop corrections to be at most at the per-mille level. We thus take 
the remaining theoretical uncertainty in T 77 to be of order 0.1%. 



6. H -> Z~/ 

Like H — > 77, the SM H — > Z7 decay partial width receives QCD corrections to the 
one-loop diagrams involving quarks. These corrections were given in Ref. 57|] including the 
full dependence on the Higgs, Z and top quark masses. They amount to less than 0.3% in 
the intermediate Higgs mass range [16( and have been neglected in HDECAY. 

In analogy with the H — > 77 process, we expect the two-loop electroweak corrections to 
T^ 7 to be at the few-percent level; we thus take an overall theory uncertainty of 4% for the 
H — > Z7 mode. 



APPENDIX B: THEORETICAL UNCERTAINTY IN THE e+e 
DUCTION CROSS SECTION 



uvH PRO- 



The full one-loop electroweak radiative corrections to Higgs production via e + e — > vvH 



have been calculated in Refs. 



3,Q, 



19]. For the experimental selection used in the 1000 
GeV centre-of-mass energy ILC studies jl5|. the processes e + e~ — > y^iH [i = e,n,T) 



29 



arising from both WW fusion and Higgsstrahlung contribute, with WW fusion dominating. 



The one- 
10% 18, 



oop electroweak corrections to these processes reduce the cross section by about 



19| in the parameter range of interest to us. 



Beyond one loop, the e + e~ — ► uPH cross section will receive additional electroweak 
and QCD corrections. We estimate the missing two-loop electroweak corrections to be 
roughly of order (one- loop) 2 . The largest electroweak corrections are due to initial-state 
radiation, which is already computed in the leading-log approximation up to order a 3 19]. 
The remaining purely weak one-loop corrections are of order 3-4% 19j , leading us to estimate 
a (one-loop) 2 uncertainty of about 0.1-0.2%. At two loops, QCD corrections will affect only 
the fermionic component of the one-loop electroweak corrections; this fermionic component 
amounts to about a 2% reduction of the total cross section in our parameter region of 
interest [18J, so we estimate the QCD correction at two loops to be of order 0.2%. At 
TeV energies one should also consider the enhanced effects that often occur in high-energy 
(E ^> mw) processes, such as logarithmic running of couplings and soft or collinear weak 
boson exchange; these effects can be at the percent level even at two loops. However, the 
t-channel weak boson fusion process of interest here is dominated by small energy transfers, 
so that we do not expect these effects to exceed the 0.1% level. For a relatively light Higgs 
mass mil ~ 120-140 GeV the Higgs self-coupling is relatively small, so that we should also 
be safe from large corrections due to Higgs self-interaction diagrams at the two-loop level. 

We thus estimate a combined theory uncertainty on the e + e _ — ► vvH cross section at 



1000 GeV centre-of-mass energy of 0.5% [44J|. We emphasize that this is just a "reasonable" 



estimate; to make this number really solid would require a serious study of the size of the 
two-loop effects mentioned above. 



APPENDIX C: UNCERTAINTIES IN m^(rn^) AND rr^(m^) 

The predictions for the Higgs partial widths into bottom or charm quark pairs are com- 
puted up to three loops in terms of the MS quark masses evaluated at the scale of the 
Higgs mass. As input, the calculation requires the bottom and charm quark MS masses at 
a starting scale such as the quark's own mass: m%(mb) and m^(m^). We summarize here 
the various ways that m^m^) and m^(m^r) have been extracted; the masses obtained are 
summarized in Tables [Villi and ITXl 
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Source m c (m c ) [GeV] 
Inclusive semileptonic B decays 

via mf and m b - m c [22] 1.224 ± 0.057 (*) 

via kinetic masses [21] 1.24 ± 0.07 
e + e~ — > hadrons 

from spectral moments [62] 1.304 ± 0.027 

more conservative error analysis [61] 1.29 ± 0.07 
Quenched lattice QCD 

from F K and D s mass [65] 1.301 ± 0.034 ± ?? 

from cc and cs state masses [66] 1.319 ± 0.028 ± ?? 
Unquenched lattice QCD 

meson spectra (prelim) [68] 1.22 ±0.09 



TABLE VIII: MS charm quark masses, m^(m^), from recent analyses (see text for details). In our 
study we use the charm mass extracted from inclusive semileptonic B decays [22] , marked with an 
asterisk (*) in the table. For quenched lattice QCD the question marks ?? denote the uncontrolled 
uncertainty due to the quenched approximation. 



Source mj(mj) [GeV] 
Inclusive semileptonic B decays 

via kinetic masses [21] 4.20 ± 0.04 (*) 
e + e~ — > hadrons 

from spectral moments [62] 4.191 ± 0.051 

more conservative error analysis [61] 4.22±0.11 
Quenched lattice QCD 

from bb and bs state masses [66] 4.33 ±0.10 ± ?? 
Unquenched lattice QCD 

T spectrum [67] 4.4 ± 0.3 

meson spectra (prelim) [68] 4.7 ±0.4 

TABLE IX: MS bottom quark masses, ra&(mb), from recent analyses (see text for details). In our 



study we use the bottom mass extracted from inclusive semileptonic B decays 21], marked with an 
asterisk (*) in the table. For quenched lattice QCD the question marks ?? denote the uncontrolled 
uncertainty due to the quenched approximation. 

Extracting the precise values of the charm and bottom quark masses from experiment 
is a challenge because of QCD effects which are starting to become strong at the bottom 
and charm mass scales. In particular, perturbative QCD does not allow one to define the 
quark pole masses with an accuracy better than Aqcd ~ 200 MeV; this is known as the 
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renormalon ambiguity. 4 Because the MS quark masses are short-distance mass definitions, 
they are in principle free of this ambiguity. 

1. Inclusive semileptonic B decays 

I — I 

The charm mass has recently been extracted [22| from a global fit to B meson decay 



spectra 



58] . The procedure combines the values of mj, — m c and ml from a fit of the decay 



spectra of B — ► X c lv and B — > to a Heavy Quark Effective Theory (HQET) expansion 



performed in Ref. 58]. The proper treatment of the renormalon cancellation was done in 



Ref. 



22] to extract the charm mass m c (m c ) = 1.224 ± 0.017 exp t ± 0.054 t h y GeV, where the 



first error is experimental and includes the uncertainty in a s and the second error represents 



a conservative combination of the theory uncertainties in the extraction [22] . While this 
theory uncertainty is not meant to be interpreted as a statistical uncertainty (it represents 
a range where the true charm mass should be located with probability much higher than 
the 67% probability of a lcr interval), we will nevertheless treat it as a la Gaussian error 
and combine it in quadrature with the experimental error for inclusion in our y 2 fits. For 
our analysis we thus take a charm mass of 

fn c {mt) = 1.224 ± 0.057 GeV. (CI) 

The quoted uncertainty amounts to a 4.7% relative uncertainty, which is dominated by the 
theoretical uncertainty. The theoretical uncertainty could be improved in the future by a 
full 0(a>1) analysis of the inclusive B decay spectra used in the extraction, to reduce the 
theoretical uncertainty in the fit to the HQET parameters. 



The charm quark mass has also been extracted in Ref. [2lJ] through a similar fit to the 
B — > X c lv and B — > X s 7 decay distributions using HQET expansions in the kinetic scheme, 
yielding kinetic masses rrib and m c which are converted to MS masses using a formula 
from Ref. 0], accurate to two loops and including part of the three-loop al corrections. 
The resulting MS charm mass is m^(rn^) = 1.24 ± 0.07 GeV. This is consistent within the 



4 A renormalon is a power-like infrared-sensitive contribution to a hard process - a relative correction of 
order Aqcd/wq in our case - as opposed to the usual logarithmic infrared-sensitive effects log(AQCD/'7iQ) 
encountered in QCD calculations. 
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uncertainties with the charm mass found in Ref. [221 ] . 

The fits of the B — ► X c lv and -B — > Xgj decay distributions to HQET expansions 
are also used to extract the bottom quark mass. The extraction can be done using any 



of a variety of different schemes for the HQET expansions [58(; however, some schemes 
are better than others in the sense that they suffer from smaller theoretical uncertainties. 
Currently the smallest uncertainty comes with the bottom mass extraction in the IS" scheme, 
m\ s = 4.68 ± 0.03 GeV {|38j], an uncertainty of only 0.64%. The bottom mass has also been 
extracted 2l|, |58j in the kinetic scheme discussed above for the charm mass extraction. As 



for the charm quark, the kinetic bottom quark mass can be translated into an MS mass 



using the formula from Ref. [59(. Following this procedure, Ref. [2l| finds for the MS mass, 

rrT b (mE) = 4.20 ± 0.04 GeV. (C2) 

We use this mass in our analysis. The relative uncertainty is 0.95%, and is dominated by 
the theoretical uncertainties that feed into the fit errors. 

2. e + e~ — > hadrons 

The inclusive cross section for e + e~ — > hadrons as a function of the centre-of-mass en- 
ergy is sensitive to the fundamental parameters of QCD, including the quark masses. In 
particular, precise perturbative QCD calculations of the first few "moments" of the cross 
section, which depend heavily on the threshold region, allow a precise determination of the 
MS charm and bottom quark masses. 

The moments P n = J R qq -(s)s~( n+1 ^ds, where R gg - = a(e + e~ — > qq + X) /a(e + e~ — > 
and q = c or b, are obtained experimentally from measurements of the production cross 
section near the cc and bb thresholds. Experimental uncertainties include gaps in the data 
above the quarkonium resonance re gion s and uncertainties in the large contributions from 



the narrow quarkonium resonances 60|. The quark masses are extracted by matching the 
experimentally determined moments to theoretical predictions. The uncertainties on the 
theory side include the unknown higher-order QCD contributions to the moments, which 
manifest both as scale dependence within a single scheme and differing results for the quark 
masses depending on how the renormalization scale is chosen as a function of the e + e~ 
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centre-of-mass energy 61] . 



A recent fit 



62j found m^(m^) = 1.304 ± 0.027 GeV and m^(m^) = 4.191 ± 0.051 GeV 



This remarkably small error estimate for m c was reconsidered in Refs. [60j, ]6l|] , with special 
attention being paid to sources of theory uncertainty, resulting in revised mass extractions 
of rn c (nQ = 1.29 ± 0.07 GeV and WTE(m) = 4.22 ± 0.11 GeV [61]. 

Very recently the first moment of the hadronic production cross section was computed at 



0(af) |63|, ]64|, significantly reducing the renormalization scale dependence. Repeating the 
analysis of Ref. [62j, Ref. (gj finds m^(m^) = 1.295±0.015 GeV and mj(mj) = 4.205±0.058 



GeV. We await a more detailed examination of the theoretical uncertainties a la Refs. 
before making use of this very small uncertainty in m c . 



60 



61] 



3. Lattice QCD 

Lattice QCD allows the extraction of the quark masses through fits of the fundamental 
QCD parameters - a s and the quark masses - to precisely-measured meson properties such 
as the kaon mass and decay constant Fr (to set the QCD scale and the light quark masses) 
and the mass of the D s meson (to provide sensitivity to the charm quark mass) (ssfl . Adding 
B mesons to the simulation allows access to the bottom quark mass 66]. The charm and 
bottom masses have been thus calculated in quenched lattice QCD, in which only virtual 
gluons are simulated and not virtual light-quark-antiquark pairs. Recent analyses find 



m(mZ) = 1.301 ± 0.034 GeV |65| or 1.319 ± 0.028 GeV [66j and rn b (m) = 4.33 ± 0.10 GeV 
66j . However, quenched lattice QCD calculations suffer from an uncontrolled theoretical 
uncertainty due to the quenched approximation. 

Unquenched lattice QCD simulations are significantly more computer-intensive and have 
only recently yielded significant results. The bottom quark mass in full (unquenched) lattice 



QCD was obtained for the first time in Ref. 67| from the T spectrum. The bare mass on the 



lattice was extracted with 1-2% uncertainty, which is mostly systematic and due to missing 
radiative corrections to the leading relativistic and discretization corrections to the T binding 
energies. The bare mass on the lattice depends on the lattice spacing, which cuts off the 
effective theory on the lattice; it can be converted to the MS mass without encountering 



infrared problems. The conversion is done at one loop yielding mi,(m&) = 4.4 ±0.3 GeV 671 ]. 



where the error is dominated by the perturbative uncertainty in the conversion from the bare 
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lattice mass. A two-loop calculation of the heavy quark self-energy to reduce the conversion 



error is underway 



6 71 ] - Reducing the error in mb(mb) to the level of the uncertainty in the 



bare lattice mass would make the lattice determination of mft(m&) competitive with that from 
inclusive B meson decay distributions discussed above. A similar (preliminary) analysis 68] 
including also the charm quark gives m^(m^) = 1.22 ± 0.09 GeV and m^(rnb) = 4.7 ± 0.4 
GeV. These masses are all in good agreement within their uncertainties with the charm and 
bottom masses found through other techniques. 
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